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A  coordinate  independent  derivation  of  the  Eulerian  and  Lagrangian  strain 
tensors  of  finite  deformation  theory  is  given  based  on  the  parallel  propaga¬ 
tor,  the  world  function,  and  the  displacement  vector  field  as  a  three-pomt 
tensor.  The  derivation  explicitly  shows  that  the  Eulerian  and  Lagrangian 
strain  tensors  are  two-point  tensors,  each  a  function  of  both  the  spatial 
and  material  coordinates.  The  Eulerian  strain  is  a  two-point  tensor  that 
transforms  as  a  second  rank  tensor  under  transformation  of  spatial  coor¬ 
dinates  and  transforms  as  a  scalar  under  transformation  of  the  material  co¬ 
ordinates.  The  Lagrangian  strain  is  a  two-point  tensor  that  transforms  as  a 
scalar  under  transformation  of  spatial  coordinates  and  transforms  as  a  sec¬ 
ond  rank  tensor  under  transformation  of  the  material  coordinates.  These 
transformation  properties  are  needed  when  transforming  the  strain  tensors 
from  one  frame  of  reference  to  another  moving  frame. 
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1.  Background 


The  U.S.  Army  is  developing  an  electromagnetic  gun  (EMG)  for  battlefield 
applications.  During  the  past  few  years,  on  a  recurring  basis.  Dr.  John  Lyons 
(ARL  Director)  and  Dr.  W.  C.  McCorkle  (Director  of  U.  S.  Army  Aviation 
and  Missile  Command)  have  requested  that  I  look  at  some  of  the  physics 
of  the  EMG.  In  the  most  recent  request,  I  was  asked  to  look  at  stresses  in 
a  rotating  cylinder.  For  an  elastic  cylinder,  this  is  a  classic  problem  that  is 
solved  in  many  texts  on  linear  elasticity  [1-6].  However,  when  these  deriva¬ 
tions  are  examined  closely,  one  finds  certain  shortcomings  [7].  Therefore,  I 
spent  some  time  looking  at  the  problem  of  stresses  in  elastic  rotating  cylin¬ 
ders,  which  resulted  in  a  manuscript  [7].  In  the  course  of  this  work  [7], 
I  had  to  clearly  imderstand  the  transformation  properties  of  the  Largan- 
gian  and  Eulerian  strain  tensors  of  finite  deformation  theory.  I  was  quite 
dissatisfied  with  the  standard  derivations  of  the  Largangian  and  Eulerian 
strain  tensors  because  these  derivations  take  either  of  two  (both  xmpalat- 
able)  approaches.  The  first  approach  uses  shifter  tensors,  which  are  often 
defined  as  inner  products  between  two  basis  vectors  at  two  different  spatial 
locations  [8,9].  In  this  approach,  basis  vectors  are  not  parallel  transported 
to  the  same  spatial  location  before  the  inner  product  is  carried  out.  This 
is  unpalatable,  even  in  Euclidean  space,  unless  one  is  using  Cartesian  co¬ 
ordinates.  The  second  approach  uses  convected  (moving)  coordinates,  and 
vectors  and  tensors  are  associated  with  a  given  coordinate  in  the  convected 
(moving)  coordinate  system,  rather  than  being  associated  with  a  point  in 
the  imderl3ting  space. 

In  the  derivation  that  I  present  below,  I  avoid  both  of  the  vmpalatable  fea¬ 
tures  mentioned  above.  I  provide  a  coordinate  independent  derivation  of 
the  Lagrangian  and  Eulerian  strain  tensors  based  on  standard  concepts  in 
differential  geometry:  the  parallel  propagator,  the  world  function,  and  the 
displacement  vector  field  as  a  three-point  tensor. 

The  derivation  that  I  present  below  is  also  useful  for  gaining  a  basic  imder- 
standing  of  the  role  of  the  xmstrained  state,  or  reference  configuration,  in 
the  definition  of  the  strain  tensors.  Having  a  firm  conceptual  grasp  of  the 
role  of  the  imstrained  state  in  the  definition  of  the  strain  tensors  is  imper¬ 
ative  for  xmderstanding  the  behavior  of  pre-stressed  materials  vmder  finite 
deformations  in  high-stress  applications,  such  as,  for  example,  in  the  elec¬ 
tromagnetic  rail  gun  [10]. 
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2.  Introduction 


The  theory  of  stresses  in  rotating  cylinders  and  disks  is  of  great  importance 
in  practical  applications  such  as  rotating  machinery,  turbines  and  genera¬ 
tors,  and  wherever  large  rotational  speeds  are  used.  In  a  previous  work  [7], 
1  gave  a  detailed  treatment  of  stresses  in  a  rotating  elastic  cylinder.  This  is  a 
classic  problem  that  is  treated  in  many  texts  on  linear  Elasticity  theory  [2-6]. 
These  treatments  linearize  the  strain  tensor  in  the  gradient  of  the  displace¬ 
ment  field,  assuming  that  these  (dimensionless)  gradients  are  small.  For 
large  angles  of  rotation,  the  quadratic  terms  (in  displacement  gradient  in 
the  definition  of  strain)  are  as  large  as  the  linear  terms,  and  consequently, 
these  quadratic  terms  carmot  be  dropped  [7].  In  Ref  [7],  I  provide  an  al¬ 
ternative  derivation  of  stresses  in  an  elastic  cylinder  that  relies  on  trans¬ 
forming  the  problem  from  an  inertial  frame  (where  Newton's  second  law 
is  valid)  to  the  co-rotating  frame  of  the  cylinder,  where  the  displacement 
gradients  are  small.  During  the  course  of  that  solution,  I  had  to  transform 
the  Lagrangian  and  Eulerian  strain  tensors  of  finite  elasticity  to  the  (non- 
inertial)  co-rotating  frame  of  reference  of  the  cylinder,  which  is  a  moving, 
accelerated  frame.  This  work  required  the  detailed  understanding  of  the 
transformation  properties  of  the  Lagrangian  and  Eulerian  strain  tensors. 

The  standard  derivation  of  these  strain  tensors  is  done  with  the  help  of 
shifter  tensors  [8,9].  Shifter  tensors  are  often  defined  in  terms  of  irmer  prod¬ 
ucts  of  basis  vectors  that  are  located  at  two  different  spatial  points  [8,9]. 
For  me,  inner  products  between  vectors  at  two  different  points  is  an  im- 
palatable  operation,  even  in  Euclidean  space.  In  order  to  compute  the  inner 
product  between  two  vectors,  the  vectors  must  first  be  parallel  transported 
to  the  same  spatial  point  (unless  we  are  using  Cartesian  coordinates,  in 
which  case  the  derivation  becomes  coordinate  specific). 

In  other  treatments,  where  shifter  tensors  are  not  employed  in  the  deriva¬ 
tion  of  strain  tensors,  convected  (moving)  coordinates  are  used,  see  for  ex¬ 
ample  [11-14].  When  using  convected  (moving)  coordinates,  the  coordi¬ 
nates  of  the  initial  imdeformed  point  and  the  deformed  point  are  the  same, 
but  the  basis  vectors  change  during  deformation.  In  derivations  of  strain 
tensors  using  convected  coordinates,  vectors  and  tensors  are  associated 
with  a  given  point  in  the  convected  (moving)  coordinate  system,  rather 
than  being  associated  with  a  point  in  the  underlying  (inertial)  space.  Ten¬ 
sors  are  absolute  geometric  objects,  and  they  should  properly  be  associated 
with  a  point  in  the  underlying  space,  and  not  a  given  coordinate,  e.g.,  in 
moving  coordinates. 

In  this  work,  I  avoid  the  unpalatable  features  of  the  strain  tensor  derivation 
mentioned  in  the  above  two  paragraphs.  I  derive  the  strain  tensors  using 
the  concept  of  absolute  tensors,  where  a  tensor  is  associated  with  a  point 
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in  the  space,  rather  than  the  coordinates  in  a  given  (moving)  coordinate 
system.  I  provide  a  coordinate  independent  derivation  of  the  Lagrangian 
and  Eulerian  strain  tensors,  where  I  keep  track  of  the  positions  of  the  basis 
vectors.  The  derivation  necessarily  uses  two-point  (and  three-point)  ten¬ 
sors  [8,9,15,16].  The  derivation  is  based  on  standard  concepts  in  differen¬ 
tial  geometry:  the  parallel  propagator  (a  two-point  tensor),  the  world  func¬ 
tion  (a  two-point  scalar),  and  the  displacement  vector  field  (a  three-point 
tensor).  This  derivation  makes  clear  the  transformation  properties  of  the 
strain  tensors  under  coordinate  transformations  from  one  frame  of  refer¬ 
ence  to  a  second  frame  that  is  moving  and  accelerated  (with  respect  to  the 
first  frame). 

The  derivation  below  of  the  Eulerian  and  Lagrangian  strain  tensors  makes 
the  transformation  properties  (e.g.,  to  a  moving  frame)  clear.  Furthermore, 
this  derivation  makes  the  role  of  the  reference  (unstrained)  configuration 
more  clear  in  the  definition  of  the  strain  tensors.  Clarifying  this  role  is  im¬ 
portant  for  appl5dng  finite  deformation  theory  to  pre-stressed  materials, 
which  are  capable  of  withstanding  higher-stress  applications,  such  as  in  ro¬ 
tating  machinery  [7,10].  Finally,  the  derivation  presented  here  allows  the 
generalization  of  the  definition  of  strain  tensors  to  the  realm  where  general 
relativity  applies  [17,18]. 
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3.  Geometric  Background 


In  Euclidean  space,  a  vector  can  be  trivially  parallel  propagated  in  the  sense 
that  after  a  round  trip  the  vector  still  points  in  the  same  direction.  In  Rie- 
mannian  space,  the  parallel  displaced  vector  is  not  equal  to  itself  after  the 
round  trip  parallel  displacement.  In  this  sense,  in  Euclidean  space  we  need 
not  distinguish  the  position  of  a  vector  because  "it  always  points  in  the 
same  direction  under  parallel  displacement,"  even  though  its  components 
may  be  different  from  point  to  point  because  the  basis  vectors,  onto  which 
we  project  the  vector,  point  in  a  different  direction  from  point  to  point. 
So,  in  Euclidean  space  the  parallel  displaced  (physical)  vector  (a  geomet¬ 
ric  object)  is  thought  to  point  in  the  same  physical  direction.  In  Riemannian 
space,  however,  the  situation  is  quite  different.  In  Riemannian  space,  a  vec¬ 
tor  that  is  parallel  displaced  will  in  general  point  in  a  different  direction. 
The  physical  test  is  to  parallel  displace  the  vector  along  a  curve  that  re¬ 
turns  to  the  starting  point.  If  there  is  non-zero  curvature,  as  measured  by 
the  Riemann  curvature  tensor,  then  upon  returning  to  its  starting  point  the 
vector  components  will  be  different  than  the  initial  vector  components  at 
the  starting  point.  So,  in  Reimannian  space,  it  is  imperative  to  specify  the 
position  of  a  vector.  In  Euclidean  space,  appropriate  to  material  deforma¬ 
tions,  I  also  keep  track  of  the  position  of  a  vector.  This  additional  care  in 
Euclidean  space,  together  with  the  transformation  properties  of  the  world 
fvmction,  leads  to  a  clearer  understanding  of  the  transformation  properties 
of  the  Lagrangean  and  Eulerian  strain  tensors,  under  transformations  from 
one  system  of  coordinates  to  another  that  is  in  relative  motion  (a  moving 
frame). 

This  section  briefly  reviews  the  fundamental  geometric  quantities  that  nat¬ 
urally  arise  in  discussion  of  deformation,  but  which  are  not  usually  dis¬ 
cussed  in  this  context.  These  quantities  are  the  world  function  (or  funda¬ 
mental  two-point  scalar  of  the  the  space),  the  parallel  propagator,  and  the 
position  vector.  This  section  will  also  serve  to  define  my  notation.  Each  of 
the  quantities  mentioned  are  examples  of  a  class  of  geometric  object  known 
as  two-point  tensors,  which  occur  naturally  m  the  discussion  of  deforma¬ 
tions.  I  have  foimd  useful  discussions  of  general  tensor  calculus  in  Synge 
and  Schild  [19]  and  Synge  [16],  and  discussions  oriented  toward  deforma¬ 
tion  theory  in  the  appendix  by  Ericksen  in  Treusdell  and  Toupin  [15],  and  in 
Narasimhan  [9],  Eringen  [8],  and  Eringen  [20].  In  particular,  discussion  of 
two-point  tensors  can  be  found  in  Synge  [16],  Ericksen  [15],  Narasimhan  [9] 
and  Eringen  [8]. 

3.1  The  World  Function 

The  world  function  was  initially  introduced  into  tensor  calculus  by 
Ruse  [21,22],  Synge  [23],  Yano  and  Muto  [24],  and  Schouten  [25].  It  was 
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further  developed  and  extensively  used  by  Synge  in  applications  to  prob¬ 
lems  dealing  with  measurement  theory  in  general  relativity  [16].  An  ap¬ 
plication  of  the  world  function  to  problems  of  navigation  and  time  transfer 
can  be  found  in  Ref.  [26],  Compared  to  the  enormous  attention  given  to  ten¬ 
sors,  the  world  function  has  been  used  very  little  by  physicists.  Yet,  when 
geometry  plays  a  central  role,  such  as  in  deformation  theory,  the  world 
hmction  is  helpful  to  clarify  and  unify  the  tmderl5dng  geometric  concepts. 
The  world  function  is  simply  one-half  the  square  of  the  distance  between 
two  points  in  the  space.  In  applications  to  relativity  and  four-dimensional 
space-time,  the  space-time  is  often  taken  as  a  general  (curved)  pseudo- 
Riemannian  space  [16].  In  applications  to  deformation  of  materials,  we  are 
concerned  with  a  Euclidean  tiuree-dimensional  space.  However,  for  under¬ 
standing  the  transformation  properties  of  displacement  vectors  and  strain 
tensors,  it  is  helpful  to  use  the  concept  of  world  function  in  a  Euclidean 
three-dimensional  space  described  by  curvilinear  coordinates  x\i  =  1,2,S, 
with  a  metric  gij,  which  in  general  is  a  function  of  position. 


Consider  two  points  in  a  general  Riemannian  space.  Pi  and  P2,  connected 
by  a  unique  geodesic  path  (a  straight  line  in  Euclidean  space)  T,  given  by 
x*(u),  i  =  1, 2, 3,  where  ui  <  u  <  U2,  and  x*(u)  are  curvilinear  coordinates 
of  the  path.  The  coordinates  of  point  Pi  =  {x\}  and  point  P2  = 
general,  a  geodesic  is  defined  by  a  class  of  special  parameters  u',u--  •,  that 
are  related  to  one  another  by  linear  transformations  u'  =  au  +  h,  where 
a  and  h  are  constants.  Here,  it  is  a  particular  parameter  from  the  class  of 
special  parameters  that  define  the  geodesic  T,  and  x*(ti)  satisfy  the  geodesic 
equations 

<Px^  dx^dx'^ 


du^ 


+  ^k-T- 
du 


du 


=  0 


(1) 


Using  Cartesian  coordinates  (rather  than  general  curvilinear  coordinates 
a;*')  in  Euclidean  space,  the  Christoffel  symbol  T*-;^  =  0,  and  the  solution  of 
equation  (1)  is  simply 


=  z?  +  {z^  -  zf)  (2) 

U2  —  Ui 

where  ui  <  u  <  U2,  i  —  1,2,3  and  the  Cartesian  coordinates  of  points 
Pi  and  P2  are  zf  and  Z2 ,  respectively.  In  a  general  Riemannin  space,  the 
world  function  between  point  Pi  and  P2  is  defined  as  the  integral  along  T 
in  arbitrary  curvilinear  coordinates  x*  by 

=  (3) 

The  value  of  the  world  fimction  has  a  simple  geometric  meaning:  It  is  one- 
half  the  distance  between  points  Pi  and  P2.  Its  value  depends  only  on  the 
eight  coordinates  of  the  points  Pi  and  P2.  The  value  of  the  world  function 
in  equation  (3)  is  independent  of  the  particular  special  parameter  u  in  the 
sense  that  vmder  a  transformation  from  one  special  parameter  u  to  another, 
u',  given  by  tt  =  au'  -f  b,  with  a;*(u)  =  x''{u{u'))  and  with  a  and  b  constant. 
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the  world  function  definition  in  equation  (3)  has  the  same  form  (with  u 
replaced  by  u'). 

The  world  function  is  the  fundamental  two-point  invariant  that  character¬ 
izes  the  space.  It  is  invariant  under  independent  transformation  of  coordi¬ 
nates  at  Pi  and  at  P2.  For  a  given  space,  the  world  function  between  points 
Pi  and  P2  has  the  same  value  independent  of  the  coordinates  used,  which 
makes  it  a  useful  coordinate  independent  quantity.  In  Euclidean  space,  us¬ 
ing  Cartesian  coordinates,  the  world  function  has  the  simple  form 

Q.{z\ ,  4)  =  i  Sij  Az"  Az^  (4) 

where  6ij  is  the  Euclidean  metric  with  only  non-zero  diagonal  components 
(-1-1,  -1-1,  -1-1),  and  Az*  =  {z^  —  zl),  i  =  1, 2, 3,  where  z\  and  z^  are  the  Carte¬ 
sian  coordinates  of  points  Pi  and  P2,  respectively.  (I  use  the  convention  that 
all  repeated  indices  are  summed,  unless  otherwise  stated.) 

The  world  function  has  a  number  of  interesting  properties,  see  Synge  [16]. 
Calculations  of  the  world  function  for  spaces  other  than  Euclidean  spaces, 
namely  four-dimensional  space-time,  can  be  found  in  Refs.  [16,26-29].  In 
what  follows,  I  restrict  myself  to  a  three-dimensional  space.  By  transform¬ 
ing  to  a  new  system  of  coordinates,  say  spherical  coordinates. 


r  COS0 

cos 

(5) 

r  cos  9 

sincf) 

(6) 

r  cos  9 

(7) 

the  world  function  in  equation  (4)  can  be  expressed  as  a  function  of  spheri¬ 
cal  coordinates  of  point  Pi  =  {ri.,6i,(f)i),  and  P2  =  (r2,  ^2,  ^2)- 

Consider  a  geodesic  given  by  equation  (1)  in  a  general  three-dimensional 
Riemannian  space.  The  covariant  derivatives  of  the  world  function  have 
two  important  properties; 


dn{Pi,P2) 

a  4 

dniPi,P2) 

dx\ 


{U2  -  Ui) 


L  Ajj 


-{U2  -  Ui) 


-LXi, 


(8) 

(9) 


where 


L  -  [2  fl(Pi,P2 )]'/"=  ds=  [ 

JP,  Ju 


r'^2 

du  du 

1/2 


du 


(10) 


is  the  length  of  the  geodesic  between  Pj  and  P2,  gij  is  the  metric  in  coordi¬ 
nates  X®,  and  Ajj  and  Aj2  are  components  of  the  imit  tangent  vectors  at  end 
points  Pi  and  P2  (assuming  non-null  geodesics  [16]): 


Ail  “ 


Ai2 


(11) 

(12) 
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where  the  relation  between  parameter  u  and  arc  length  s  is  given  by  equa¬ 
tion  (10).  In  equations  (8)  and  (9),  the  covariant  partial  derivatives  with 
respect  to  x\  and  are  done  with  respect  to  the  coordinates  of  points  Pi 
and  P2.  See  figure  1. 

For  the  special  case  of  interest  in  deformation  of  materials,  the  space  is 
three-dimensional  Euclidean,  the  geodesic  is  a  straight  line,  and  the  vec¬ 
tors  Ajj  and  Ajj  are  colinear,  although  I  still  consider  them  as  existing  at 
distinct  points.  Using  Cartesian  coordinates  and  the  explicit  form  of  the 
world  fvmction  in  equation  (4)  and  equations  (8)  and  (9)  take  the  form 

=  ^ij  (4  -  4)  (13) 

ni,=-Sij  (4-4)  =  -L\i,  (14) 

where  I  used  Synge's  short-hand  notation  for  the  components  of  the  co¬ 
variant  partial  derivatives  by  putting  subscripts  on  the  indices  to  indicate 
which  coordinates  were  differentiated.  This  short-hand  notation  is  particu¬ 
larly  convenient  to  show  the  transformation  properties  of  the  world  fimc- 
tion  and  to  indicate  the  spatial  location  of  vectors  and  tensors.  For  example, 
the  quantity  Qi^  transforms  as  a  vector  imder  coordinate  transformations  at 
P2  and  as  a  scalar  under  coordinate  transformations  at  point  Pi.  The  quan¬ 
tity  Ajj  is  a  vector  located  at  point  P2.  Note  that  by  virtue  of  their  definitions 
in  the  left  side  of  equations  (13)  and  (14),  the  right  sides  are  two-point  ten¬ 
sors,  whose  components  are  functions  of  coordinates  at  point  Pi  and  P2.  For 
example,  the  right  side  of  equation  (14)  is  a  product  of  a  two-point  scalar 
L,  and  a  one  point  vector,  Aij  at  Pi. 


af2(Pi,P2) 

dz\ 

dn{Pi,P2) 

dz\ 


Figure  1.  A  geodesic 
path  is  shown  in  three 
dimensions  with 
tangent  unit  vector  at 
the  ends. 
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3.2  Parallel  Propagator 

Given  a  vector  with  components  at  point  Pi,  the  vector  is  said  to  be 
parallel  propagated  from  Pi  to  P2  along  a  geodesic  curve  C  specified  by 
x'(u),  ui  <  u  <  U2,  where  Pi  =  {a;®(ui)}  and  P2  =  {x^{u2)},  when  its 
covariant  derivative  is  zero  along  this  curve; 


dv'‘ 

du 


•  dx^ 
du 


=  0 


(15) 


Equation  (15)  is  a  mapping:  Given  the  components  of  a  vector,  v*'  at  point 
Pi,  we  obtain  the  components  of  the  parallel  transported  vector  at  point 
P2  by  solving  equation  (15).  It  is  convenient  to  define  a  two-point  tensor, 
called  the  parallel  propagator  [16],  which  gives  the  components  of  a 
vector  under  parallel  translation  of  the  vector  from  point  Pi  to  point  P2. 
Given  a  vector  with  components  v*'  at  point  Pi,  the  propagator  relates 
the  components  of  this  vector  at  Pi  to  the  components  of  this  same  vec¬ 
tor  after  parallel  translation  to  point  P2 

v*'  (16) 


In  a  general  Riemannian  space,  the  components  of  the  vector  at  point  P2 
depend  on  the  path  of  parallel  translation  from  Pi  to  P2,  in  the  sense  that 
the  path  must  be  a  geodesic  by  the  definition  of  the  parallel  propagator. 
However,  in  Euclidean  space  these  components  are  completely  path  inde¬ 
pendent;  the  components  depend  only  on  the  end  points  Pi  and  P2. 

A  vector  is  considered  as  a  geometric  object,  which  means  that  it  is  inde¬ 
pendent  of  coordinate  system.  In  a  Riemannian  space,  imder  the  operation 
of  parallel  propagation  a  vector  changes  in  such  a  way  that  its  magnitude 
stays  the  same  but  its  absolute  direction  can  change  because  of  the  curva¬ 
ture  of  the  space  [30].  The  direction  of  the  parallel  propagated  vector  is,  of 
course,  referred  to  the  local  basis  vectors.  That  the  vector  direction  changes 
under  parallel  translation  can  be  understood  by  taking  a  vector  at  point  P 
and  parallel  translating  it  over  a  curve  that  returns  to  point  P.  When  com¬ 
pared  at  point  P,  the  components  of  the  initial  vector  and  the  round-trip- 
parallel-transported  vector  will  (in  general)  be  different.  It  is  in  this  sense 
that  a  vector  changes  its  direction  under  parallel  transport. 

As  mentioned  above,  the  change  in  the  vector  that  results  under  parallel 
transport  depends  on  the  path  of  parallel  propagation  (a  geodesic).  Two 
vectors  that  are  parallel  propagated  along  the  same  path  will  maintain  the 
angle  between  them  along  the  path. 

In  a  Euclidean  space,  a  vector  (the  geometric  object)  is  considered  to  be 
unchanged  when  parallel  propagated.  The  only  thing  that  happens  is  that 
the  components  of  the  vector  on  the  local  basis  must  change  according  to 
what  is  required  to  keep  the  vector  "pointing  in  the  same  direction." 

In  Euclidean  space,  the  parallel  propagator  in  Cartesian  coordinates  is 
trivial — its  components  are  just  the  components  of  a  delta  function.  The 
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components  of  a  vector  at  point  Pi  are  related  to  the  components  of  the 
same  vector  parallel  translated  to  point  P2  by  the  propagator  (whose  com¬ 
ponents  are  given  in  a  Cartesian  coordinate  basis): 

5^2  1+1  *  =  i  (17) 

1  0  i^j 

Equation  (17)  agrees  with  our  notion  from  elementary  geometry  that  in 
Cartesian  coordinates  the  vector  components  are  constant  tmder  parallel 
propagation.  However,  using  the  parallel  propagator  in  Cartesian  coordi¬ 
nates,  we  can,  for  example,  compute  the  propagator  in  curvilinear  co¬ 
ordinates  =  (r,  6, 4>)  given  in  equation  (5)-(7),  by  the  two-point  tensor 
transformation  rule 


.  dx\P2)  dz-{P{) 

y  31  Qzm  dxi 


(18) 


The  parallel  propagator  5*  is  a  two-point  tensor  because  it  transforms  as 
a  vector  xmder  coordinate  transformation  at  point  Pi  and  xmder  coordinate 
transformation  at  point  P2. 

In  Cartesian  coordinates,  when  the  points  are  made  to  coincide,  P2  ^  Pi, 
the  propagator  reduces  to  a  Kronecker  delta  at  point  Pi;  5'^^^  (Pi). 

In  general  curvilinear  coordinates,  when  the  points  Pi  and  P2  coincide,  the 
parallel  propagator  reduces  to  the  mixed  components  of  the  metric  tensor 

lim  g\-yg^j{Pi)  (metric  at  Pi)  (19) 

The  mixed  components  of  the  metric  tensor  at  Pi,  g^j{Pi)  =  g^^{Pi)  QkjiPi)/ 
are  a  Kronecker  delta — a  unit  tensor  whose  components  are  the  same  in  all 
systems  of  coordinates.  Indices  can  be  lowered  on  two-point  tensors  using 
the  appropriate  metric.  For  example,  the  index  i  of  the  propagator  g^^^^  can 
be  lowered  by  using  the  metric  tensor  at  point  P2: 

9k2ji  =  9kiiP2)  9\  (20) 

When  the  points  are  made  to  coincide,  P2  — >  Pi,  the  covariant  components 
of  the  propagator  become  ttie  covariant  components  of  the  metric  tensor  at 
Pi,  9k2ji  9kj{Pi),  where  9kj{Pi)  is  the  metric  at  Pi. 

The  covariant  derivatives  of  the  world  function  fl(Pi,  P2)  between  points 
Pi  and  P2  are  related  to  the  parallel  propagator  by  [16] 


^hji  ~  9iiji 

(metric  at  Pi) 

(21) 

^hj2  “  ^j2h 

=  -9iij2  =  -9j2h  (paraUel  propagator) 

(22) 

^i2j2  —  9i2j2 

(metric  at  P2) 

(23) 

Other  useful  properties  of  the  parallel  propagator  are  discussed  by  Synge  [16]. 
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3.3  Position  Vector 


The  position  vector  occupies  a  central  role  in  deformation  theory.  For  this 
reason,  I  discuss  it  in  detail  below.  In  elementary  geometry,  a  point  P  can 
be  identified  by  its  position  vector  r,  which  can  be  specified  in  Euclidean- 
space  Cartesian  coordinates  as 


r  =  (24) 

where  z”  are  the  Cartesian  components  of  the  vector  r  and  also  the  Carte¬ 
sian  coordinates  of  the  point  P.  In  terms  of  general  coordinate  basis  vectors 
e„  =  d/dx"  associated  with  the  curvilinear  coordinates  x",  the  vector  r  is 
given  by 

r  =  z^A\{P)emiP)  =  CemiP)  (25) 

The  position  vector  is  a  geometric  object  at  point  P.  Among  all  basis  vec¬ 

tors,  the  Cartesian  basis  vectors  in  are  unique  in  that  they  are  usually  not 
associated  with  a  particular  spatial  point.  However,  when  we  express  these 
Cartesian  basis  vectors  in  terms  of  curvilinear  basis  vectors  Bn,  then  we 
must  imagine  that  these  basis  vectors  exist  at  a  particular  point  P.  Hence, 
the  transformation  between  the  Cartesian  basis  vector  and  curvilinear 
coordinate  basis  vectors  at  point  P  associated  with  coordinates  x*  is 
given  by 

iniP)  =  A\iP)em{P)  (26) 

where  the  matrix  A™(P)  depends  on  the  coordinates  of  point  P: 

■4”  (^’)  =  ^(P)  (27) 

In  Cartesian  coordinates,  the  components  of  the  vector  r  are  simply  the 
Cartesian  coordinates  z"  of  point  P.  The  three  numbers  (z\  z^,z^)  trans¬ 
form  as  the  components  of  a  vector  vmder  orthogonal  coordinate  transfor¬ 
mations.  Note  that  in  curvilinear  coordinates,  the  components  of  the  po¬ 
sition  vector,  are  not  the  curvilinear  coordinates  of  point  P.  Also,  note 
that  the  position  vector  r  of  point  P  has  a  magnitude  equal  to  the  Euclidean 
length  from  the  origin  of  coordinates,  say  point  O,  to  point  P.  The  position 
vector  of  point  P  is  a  geometric  object  at  point  P,  however,  it  also  depends 
on  the  point  O.  This  dependence  on  point  O  is  coordinate  independent. 
Therefore,  the  position  vector  of  point  P  is  a  two-point  tensor;  it  depends 
on  point  P  and  on  point  O.  The  transformation  properties  of  the  position 
vector  are  that  of  a  scalar  when  a  change  of  coordinates  is  made  at  point  O 
and  the  transformation  is  that  of  a  vector  when  coordinates  at  point  P  are 
changed. 

In  a  Riemannian  (a  generalization  of  Euclidean  space)  space,  the  compo¬ 
nents  of  the  position  vector  ri{P)  at  point  P  can  be  defined  in  terms  of  the 
covariant  derivative  of  the  world  function 

np  =  =  ^ip{0,  P)  =  [2  f^(0,  P)]'/2  h{P)  (28) 

O  Zp 
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where  r^p  is  a  xmit  vector  at  point  P  and  [2  Cl{0,  is  the  length  of  the 
geodesic  from  point  O  to  point  P.  For  Euclidean  space,  [20(0,  is 

die  length  of  the  straight  line  OP.  Equation  (28)  shows  explicitly  that  the 
position  vector,  rip,  is  a  two-point  tensor. 

3.4  Displacement  Vector 

Consider  an  elastic  body  that  undergoes  a  finite  deformation  in  time.  The 
deformation  can  be  specified  by  a  flow  function  or  displacement  mapping 
fimction 

z'^  =  z'^{Z^,t)  (29) 

where  the  coordinates  (here  taken  to  be  Cartesian)  of  a  particle  at  point 
Q  at  time  t  are  given  in  terms  of  the  particle's  coordinates  of  point  P 
in  some  reference  state  (configuration)  at  time  t  =  to,  so  that  z^{Z'^,to)  = 
Z^.  We  assume  the  deformation  mapping  function  has  an  inverse,  which  is 
quoted  here  for  later  reference 

Z'^  =  Z"^(2^^)  (30) 

We  assume  that  both  the  coordinates  z^  and  Z'^  refer  to  the  same  Cartesian 
coordinate  system.  * 

In  deformation  theory,  the  initial  position  of  the  particle  in  the  medium  at 
point  P  is  specified  by  a  vector 

R(P)  =  Z''ifc(P)  (31) 

and  the  final  position  is  specified  by  a  position  vector 

r{Q)  =  z'^UQ)  (32) 

where  the  quantities  ifc  (P)  and  i^  (Q)  are  the  Cartesian  basis  vectors  at  point 
P  and  point  Q,  respectively.  Conventionally,  the  deformation  of  a  medium 
is  described  by  specifying  the  displacement  "vector  field,"  which  is  defined 
as  a  difference  of  these  two  position  vectors.  However,  the  basis  vectors 
ife  (P)  and  (Q)  are  at  different  points  in  the  space.  Since  vectors  can  be  sub¬ 
tracted  only  if  they  are  at  the  same  point,  we  must  parallel  translate  ife  (P)  to 
point  Q,  or,  parallel  translate  ife(Q)  to  point  P.  Depending  on  which  map¬ 
ping  we  choose,  we  arrive  at  the  Eulerian  or  the  Lagrangian  displacement 
vector. 

First,  we  parallel  translate  vector  R(P)  to  point  Q,  and  then  subtract  the 
vectors  at  point  Q  (see  fig.  2).  This  procedure  defines  the  components  of 
the  Eulerian  displacement  vector  at  point  Q: 

u(g)  =  r(Q)-R(Q)  (33) 

*I  use  the  notation  z'°  and  for  the  Cartesian  coordinates  and  also  for  the  deformation 
mapping  functions.  There  is  no  chance  of  confusing  these  two  since  the  context  makes  it 
clear  which  one  is  used  in  a  given  instant. 
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Figure  2.  The  initial  and 
final  position  vectors, 

R(P)  and  r{Q), 
respectively,  are  shown 
as  well  as  their 
respective  parallel 
translated  vectors,  R(P) 
and  r((5).  Also,  shown 
by  a  solid  curve  is  the 
actual  displacement 
path  of  a  representative 
particle  of  the  medium, 
labeled  by 

=z^{Z^,t).The 
dashed  straight  line  is 
the  line  (geodesic) 
connecting  the  initial 
and  final  particle 
positions.  The  Eulerian 
displacement  vector  is 
u  =  r(Q)  -  R(Q),  which 
is  the  difference  of  two 
position  vectors  at  point 
Q.  The  Lagrangean 
displacement  vector, 

U  -  r(P)  -  R(P),  is  the 
difference  of  two 
position  vectors  at  point 
P. 

This  Eulerian  displacement  vector  in  equation  (33)  is  often  called  the  dis¬ 
placement  vector  in  the  spatial  representation  [8,9,31].  Alternatively,  we  can 
parallel  translate  the  vector  r((5)  to  point  P,  and  then  subtract  the  vectors 
at  point  P.  This  procedure  defines  the  Lagrangian  displacement  vector  at 
point  P: 

U(P)  =  r(P)  -  R(P)  (34) 

This  Lagrangian  displacement  vector  in  equation  (34)  is  often  called  the 
displacement  vector  in  the  material  representation  [8,9,31].  Equations  (33) 
and  (34)  show  that  these  two  vectors  are  actually  referred  to  basis  vectors  at 
different  points.  In  fact,  the  two  vectors  u((5)  and  U(P)  are  related  by  par¬ 
allel  translation.  In  a  Euclidean  space,  these  vectors  are  the  same  geometric 
objects  but  they  are  expressed  in  terms  of  basis  vectors  located  at  different 
positions  P  and  Q  (see  below). 

In  order  to  further  clarify  the  transformation  properties  of  these  two  dis¬ 
placement  vectors,  we  use  the  position  vector  as  discussed  in  the  previous 
section.  Consider  the  deformation  mapping  function  in  curvilinear  coordi¬ 
nates,  This  function  specifies  the  coordinates  (point  Q)  of  a 

particle  at  current  time  t  in  terms  of  the  coordinates  (point  P)  of  the 
particle  in  the  reference  configuration  at  time  t  =  to,  so  that 

x^{X^,to)  =  X^  (35) 
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In  addition,  tiiere  exists  a  straight  line  (a  geodesic)  F  connecting  the  points 
P  and  Q. 

The  covariant  components  of  the  position  vector  of  point  P,  R(P)  = 
BniP),  in  curvilinear  coordinates  x®  are  given  by  (see  equation  (28)) 

=  P)  =  [2  n{0,  P)f^  Rip  (36) 

OXp 

where  Rip  are  the  components  of  the  unit  vector  at  point  P  tangent  to  F 
that  cormects  point  P  and  Q.  Similcurly,  the  covariant  components  of  vector 
r((5)  =  r”  Gn{Q)  in  curvilinear  coordinates  x*  are  given  by 

=  25^  =  n,„(O.0)  =  |2fJ(0,Q)|''"  (37) 

OXq 

where  are  the  components  of  the  unit  vector  at  point  Q  tangent  to  F. 
From  equations  (36)  and  (37)  it  is  clear  that  both  Rip  and  are  two-point 
tensor  objects.  The  quantity  Rip  depends  on  point  O  and  P  and  transforms 
as  a  vector  imder  coordinate  transformations  at  P  and  as  a  scalar  tmder 
coordinate  transformations  at  point  O.  The  quantity  rj^  transforms  as  a 
vector  imder  coordinate  transformations  at  point  Q  and  as  a  scalar  under 
coordinate  transformations  at  point  O. 

The  components  of  the  Eulerian  displacement  vector  in  equation  (33)  (at 
point  Q)  are  defined  in  terms  of  the  components  of  R^^  parallel  translated 
to  point  Q: 

(38) 

=  9%  n^'^{0,P)  (39) 

where  (O,  P)  =  (P)  fifcp  {O,  P),  and  gj'^{P)  is  the  metric  at  point  P  in 

coordinates  x*.  The  contravariant  components  of  the  Eulerian  displacement 
vector  in  equation  (33)  are  given  by 

(O,  Q)  -  g\  (O,  P)  (40) 

Similarly,  the  components  of  the  Lagrangian  displacement  vector  are  given 

=  g\  (O,  Q)  -  (O,  P)  (41) 

The  vectors  whose  components  are  and  are  related  by  parallel 
transport  along  the  geodesic  F  connecting  P  and  Q  (not  along  the  parti¬ 
cle  displacement  line  given  by  equation  (29)).  Transporting  to  point  Q 


=  9\ 

(42) 

=  9%  [9^lQ^'^^{o,Q)-n^no,P)] 

(43) 

=  s\{Q)n’^Ho,Q)  - 9%n^no,P) 

(44) 

= 

(45) 
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where  in  the  transition  from  equation  (43)  to  (44)  I  have  used  the  identity 
satisfied  by  the  parallel  propagator: 


i\(Q)  =  s‘%  9^1,  (46) 

where  S\{Q)  is  a  unit  tensor  (delta  function)  at  point  Q.  Equation  (46)  is 
the  statement  that  parallel  propagation  of  a  vector  has  an  inverse,  so  the 
result  that  U^{P)  and  u\Q)  are  related  by  parallel  transport  is  true  in  both 
Euclidean  and  Reimannian  spaces. 
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4.  Strain  Tensors 


At  time  t  —  to,  consider  two  particles  in  the  medium  that  are  at  posi¬ 
tions  Pi  and  P2,  respectively,  and  are  separated  by  a  finite  distance  AS  = 
[2Q{Pi,P2)]^^^.  At  a  later  time  t  >  to  these  particles  have  moved  to  new 
positions  Qi  and  Q2  and  are  separated  by  a  distance  As  =  [2 
(see  fig.  3). 

As  a  measure  of  strain,  I  take  the  4-pomt  scalar 

^{Pi,P2-,Qi,Q2)  =  (As)2  -  {ASf  =  2n{Qi,Q2)  -  2f2(Pi,P2)  (47) 

Note  that  4’(Pi,  P2;  Qi,  Q2)  depends  on  foiur  points  Pi,  P2,  Qi,  and  Q2,  and 
by  virtue  of  its  definition  in  terms  of  the  world  ftmction,  ^  is  a  true  four- 
point  scalar  invariant  under  separate  coordinate  transformations  at  each  of 
these  four  points.  In  Cartesian  coordinates,  equation  (47)  becomes 

zi,  zi)  =  Sij  {zi  -  4)  (4  -  4)  -  Sij  {Zi  -  Zi)  {Z{  -  zi)  (48) 

where  z\  —  z\Z'^,t)  and  =  z\Z^,t)  and  they  are  related  to  the  refer¬ 
ence  configuration  at  t  —  to  by 

z\Zr,to)=-Zi  (49) 


Figure  3.  The  position 
Pi  and  P2  of  two 
particles  is  shown  in  the 
reference  configuration 
at  t  =  to,  and  the  ^  ^  ^0 

positions  <5i  and  Q2  of 
the  same  two  particles 
is  shown  at  later  time 
t  >  to-  The  path  of  the 
particles  is  shown  in 
solid  lines  and  their 
displacement  is  shown 
in  dashed  lines. 
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and 


z\Z^,to)  =  Z^^  (50) 

and  z’'{Z^,t)  is  the  deformation  mapping  function  in  Cartesian  coordi¬ 
nates,  given  in  equation  (29).  If  the  particle  positions  Pi  and  P2  are  con¬ 
sidered  as  separated  by  an  infinitesimal  distance,  then,  by  assuming  conti¬ 
nuity  in  the  medium  and  a  finite  time  t  —  to,  the  points  Qi  and  Q2  are  also 
infinitesimally  separated.  However,  because  t  — to  is  finite,  the  distance  be¬ 
tween  Pi  and  Qi  is  finite  (not  infinitesimal).  Expanding  equation  (50)  about 
the  initial  position  of  the  first  particle 

4  =  z^Zl  t)  +  t)  {Zi  -zi)  +  --.  (51) 

using  zl  =  z^{Zi,  t),  leads  to  the  relation  between  spatial  (Eulerian)  coordi¬ 
nates  z*  and  material  (Lagrangian)  coordinates  Z^ 

Az^  =  §^iZr,t)AZ^  +  ---  (52) 

where  Az®  =  z^—zl  and  AZ‘  =  Z^—Z}.  Using  equation  (52)  in  equation  (48) 
leads  to  the  measure  of  strain  in  Cartesian  coordinates 

(A42  _  (A5)2  =  ^5,,.  ||r  -  <^mn)  AZ^  AZ"  -h  •  ■  •  (53) 

=  2PmnAZ’"AZ’^  +  ---  (54) 

where  the  quantity  in  parenthesis  is  the  Lagrangean  strain  tensor,  Bmn-  The 
higher  order  terms  in  AZ’^  are  small  because  we  can  always  choose  the 
two  initial  points  Pi  and  P2  arbitrarily  close  together.  From  equations  (53) 
and  (54),  it  is  clear  that  the  Lagrangean  strain  tensor  is  a  two-point  tensor, 
depending  on  initial  point  P  (in  the  reference  configuration  at  t  =  to)  and 
point  Q  (at  time  t).  Note  that  in  equation  (54)  there  is  no  restriction  to  short 
times  t  —  to,  since  we  can  always  choose  AZ"  sufficiently  small. 

The  Eullerian  tensor  can  be  obtained  from  equation  (54)  by  using  the  fact 
that  the  flow  function  in  equation  (29)  has  an  inverse.  Using 

fiyn 

AZ-^^{z\t)AP  (55) 

We  get  the  measure  of  strain  in  terms  of  the  Eulerian  strain  tensor  e-mn- 

(r)  7^  f)  7^  \ 

Az-Az”  +  ...  (56) 

=  2emnAz^Az^  +  ---  (57) 

4.1  Strain  Tensor  Derivation  in  Curvilinear  Coordinates 

I  return  to  the  definition  of  the  measure  of  strain  given  in  equation  (47). 
In  the  reference  configuration  at  t  =  to,  consider  two  particles  at  points  Pi 
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and  P2  with  curvilinear  coordinates  Xi  and  X2.  At  a  later  time  t,  these  two 
particles  are  at  positions  Qi  and  Q2  with  curvilinear  coordinates  xi  and 
X2,  respectively.  Consider  the  first  term  on  the  right  side  of  equation  (47), 
which  is  an  integral  along  a  straight  line  QiQ2- 

1  r'U’2 

<52)  =  kw2-  wi)  /  gijU^  W  dw  (58) 

where  17®  =  dx^{w) /dw  and  where  the  geodesic  (straight  line)  is  parametrized 
by  x\w),  with  wi  <  w  <  W2  and  the  end  points  are  given  by  xi  =  x'^iwi) 
and  X2  =  x^{w2)  (see  fig.  3).  The  fimction  x^{w)  is  a  solution  of  the  geodesic 
equation  (1).  In  the  case  of  Eulidean  space,  and  assuming  points  Pi  and  P2 
are  arbitrarily  close,  the  geodesic  in  equation  (58)  is  a  straight  line  given  by 

x‘^{w)  =  x\-\-  — — ^  {x\  ~  ^1) 

W2  —  WI 

The  flow  function  in  equation  (29)  maps  the  points  Pi  and  P2  into  the 
points  Qi  and  Q2.  The  points  Qi  and  Q2  depend  on  time  t.  WiHx  reason¬ 
able  continuity  assumptions,  and  the  straight  line  given  in  equation  (59) 
with  =  k{x^  -  x{)  =  k  Ax^  and  k  =  {w2-  wi)~'^,  the  world  fimction  in 
equation  (58)  can  be  approximated  by 

1  r'^2 

^{Qi,Q2)  =  7:  {w2  -  wi)  gijiQi)  kAx"  kAx^  /  dw  (60) 

Z  Jw\ 

=  \gij{Qi)Ax^Ax^  (61) 

Similarly,  the  second  term  on  the  right  side  of  equation  (47)  can  be  approx¬ 
imated  as 

fl(Pi,  P2)  =  ^  gaiPi)  AX®  AX^  (62) 

where  the  coordinates  X®  are  the  undeformed  ones  and  AX®  =  X2  —  XJ. 
The  measure  of  strain  in  equation  (47)  is  then 

{Asf  -  {ASf  =  gijiQ)  Ax®  Ax^  -  gij{P)  AX®  AX^’  (63) 

or  using  the  flow  function  in  equation  (29), 

(A»)"-(AS)^=  ^9«(0)^(P,Q)||i(P,(3)-m(P))  AX‘AX‘  (64) 

Note  that  x®  and  X®  refer  to  the  same  system  of  coordinates.  I  have  dropped 
tile  subscripts  on  Q  and  P  since  Qi  and  <32/  and  Pi  and  P2  are  infinitesi¬ 
mally  close,  respectively.  The  quantity  in  parenthesis  on  the  right  side  of 
equation  (64)  is  twice  the  Lagrangian  strain  tensor: 

2E,,,,=9«{(3)^(P,Q)|^(P,Q) -Su{P)  (65) 

From  equation  (65),  it  is  clear  that  the  Lagrangian  strain  tensor  is  a  two- 
point  tensor.  Under  transformation  of  coordinates  at  point  P,  Ekpip  is  a 
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second  rank  tensor,  while  under  transformation  of  coordinates  at  point  Q, 
it  is  a  scalar.  The  deformation  gradient,  dx^ /dX^,  is  itself  a  two-point  tensor, 
as  can  be  seen  by  its  transformation  property  when  coordinates  at  P  and  Q 
are  changed. 

It  is  interesting  to  note  that  the  Lagrangian  strain  tensor  Eki  is  convention¬ 
ally  thought  to  be  a  fvmction  of  material  coordinates,  X\  which  coincide 
with  the  point  P  (in  the  reference  configuration)  [8,9].  The  tensor  Ekpip  can 
be  taken  to  be  a  function  of  only  the  material  coordinates  by  using  the  flow 
mapping  function  in  equation  (29),  which  provides  a  mapping  between  all 
points  P  and  their  images,  points  Q,  imder  the  deformation.  However,  I  do 
not  pursue  this  interpretation  below. 

The  first  term  in  equation  (65)  is  the  Green  deformation  tensor: 

Cki  =  gijiQ)  ^(P,  Q)  |^(P,  Q)  (66) 

The  point  P  is  in  the  reference  configuration  at  time  t  —  to  and  the  point 
Q  is  in  the  deformed  state  at  time  t.  In  equation  (66),  the  Green  deforma¬ 
tion  tensor  is  naturally  a  second  rank  tensor  with  respect  to  material  coor¬ 
dinate  transformations  (point  P)  and  it  is  a  scalar  with  respect  to  spatial 
coordinate  transformations  (point  Q).  However,  the  Green  tensor  is  con- 
ventionally  taken  [8,9]  as  a  function  of  material  coordinates  by  using  the 
flow  mapping  function  in  equation  (29). 

Returning  to  equation  (64),  and  using  the  flow  mapping  function  in  equa¬ 
tion  (29),  we  can  obtain 

9ijiQ)-gki{P)-^{P,Q)-^{P,Q)j  (67) 

where  the  Eulerian  strain  tensor  is  given  by 

BX^  BX^ 

=  gij{Q)  -  gki{P)  -^{P)  ^(P)  (68) 

The  Eulerian  strain  tensor  et^j^  is  a  two-point  tensor  that  is  second  rank 
with  respect  to  spatial  coordinates  at  point  Q  and  is  a  scalar  with  respect 
to  material  coordinates  at  point  P.  Once  again,  however,  using  the  flow 
mapping  function  in  equation  (29),  ei^j^  can  be  imagined  to  depend  on  the 
material  coordinates  x®  (point  Q)  of  the  deformed  state. 


5.  Strain  Tensors  in  Terms  of  the  Displacement  Field 


The  Eulerian  and  Lagrangian  strain  tensors  can  be  expressed  in  terms  of 
the  displacement  fields  and  in  equations  (40)  and  (41).  In  terms  of 
covariant  components,  the  displacement  field  vector  at  Q  is  given  by 

(O,  Q)  -  (O,  P)  (69) 

The  covariant  derivative  of  with  respect  to  point  Q  (coordinates  x^)  is 

dX^ 

~  ^jplpiO,  P)  (70) 

where  I  used  the  chain  rule  for  covariant  differentiation  since  Q{0,  P)  = 
fl{0,P{X))  and  the  material  coordinates  X’'  =  X‘(x^)  are  functions  of 
the  spatial  coordinates  x'^.  It  is  clear  that  the  chain  rule  must  be  used  in 
equation  (70)  by  considering  Cartesian  coordinates.  In  equation  (70),  we 
also  used  the  fact  that  in  Euclidean  space  the  parallel  propagator  is  a  con¬ 
stant  imder  covariant  differentiation.  The  notation  Qi^kgiO,  Q)  means  the 
{i,  k)  component  of  the  second  covariant  derivative  of  ti\e  world  fimction  at 
point  Q.  In  Euclidean  space,  these  second  covariant  derivatives  are  simply 
related  to  the  parallel  propagator  (see  eqs  (21)-(23)). 

Using  equation  (23),  the  covariant  derivative  of  the  displacement  field  in 
equation  (70)  becomes 

dX^ 

~  9ik(Q)  —  9iQ^  9jl{P) (71) 

=  9ik{Q)  —  9iQlp  (72) 

where  in  the  last  line  the  metric  at  P,  gji  (P),  was  used  to  lower  the  index  on 
the  propagator.  Now  we  multiply  equation  (72)  by  the  parallel  propagator 
gmpiq^  sum  on  iq,  and  solve  for  the  inverse  displacement  gradient 

^  (73) 

From  equation  (73),  it  is  clear  that  in  Euclidean  space,  the  deformation  gra¬ 
dient  dX'^/dx^  is  simply  related  to  the  covariant  derivative  of  the  displace¬ 
ment  field,  Uig-kq-  Note  however,  that  in  a  Riemannian  space,  for  finite  de¬ 
formations,  it  is  generally  not  possible  to  solve  for  the  deformation  gradi¬ 
ent  [18].  Now,  using  equation  (73)  in  equation  (67),  we  find  the  expression 
relating  the  two-point  Eulerian  strain  tensor  a^jg  =  Cigjg  (P,  Q)  to  the  co¬ 
variant  derivatives  of  the  three-point  displacement  field 

(As)^  —  (A.?)^  =  '^iq\jq '^jq-,iq  ~  9  iQ) '^kq-,iq '^lq-,jq  Ax  Ax-^  (74) 
=  2eigjg  Ax'Ax^  (75) 
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where 


(76) 


^iqjQ  ~ 


1  r 

2 


9  {Q)  '^kn^in 


q;iq  '^Iq]Jq 


Equation  (76)  explicitly  shows  that  the  Eulerian  strain  tensor  is  a  func¬ 
tion  of  two  points,  material  coordinates  at  point  P  and  spatial  coordinates 
at  point  Q.  Note  that  is  not  a  function  of  point  O,  since  by  equa¬ 
tion  (71)  the  covariant  derivative  does  not  depend  on  point  O.  From 
equation  (76)  it  is  also  clear  that  transforms  as  a  second  rank  tensor 
under  spatial  coordinate  transformations  and  that  it  transforms  as  a  scalar 
under  material  coordinate  transformations. 


An  analogous  relation  can  be  obtained  for  the  Lagrangian  strain  tensor  by 
considering  the  covariant  derivative  of  the  displacement  field 


Q)-^  -  ^irkAO.P)  (77) 

Using  the  relations  between  the  second  covariant  derivatives  of  the  world 
function  and  propagator  in  equations  (21)-(23),  and  solving  for  the  dis¬ 
placement  gradient  we  get 

^  =  5 T  -  a"'"’'-  V.,:,,  (78) 

Inserting  the  displacement  gradient  in  equation  (78)  into  equation  (64)  gives 
an  expression  for  the  Lagrangian  strain  tensor  Empnp  =  Ernpnp{P’)  Q)^ 


(As)2  -  (A5)"  =  [Ump-,np  +  Unp.mp  +  9^^  {P)  U^p■,mp  Ujp,np\  AA"^  AX" 
=  2E„,pnp  AX*AX^ 


(79) 

(80) 


where 


E, 


771  pTl  p 


Ump-,np  +  Unp;mp  "i  5  {P)  Uip',mp  Pjp;np 


(81) 


Equation  (79)  shows  that  the  Lagrangian  strain  tensor  Empnp  is  a  function 
of  the  material  coordinates  at  point  P  and  spatial  coordinates  at  point  Q. 
The  Lagrangian  strain  tensor  transforms  as  a  scalar  under  spatial  coordi¬ 
nate  transformations  at  point  Q  and  as  a  second  rank  tensor  with  respect  to 
material  coordinate  transformations  at  point  P.  Note  that  there  are  sign  dif¬ 
ferences  in  equations  (76)  and  (81).  Finally,  comparing  equations  (74)  and 
(79),  we  have  the  well-known  relation  between  the  two  strain  tensors 


17  _  _  dx^  dx^ 

J^mpnp  ^iqjq  g^m  g^n 

Equation  (82)  provides  a  complicated  relation  between  the  two  two-point 
strain  tensors.  While  the  displacement  fields  u*®  and  are  related  to  each 
other  by  parallel  transport  (see  eqs  (40)  and  (41)),  the  strain  tensors  Empnp 
and  eiqjQ  are  related  by  two-point  deformation  gradient  tensors,  dx^jdX'^, 
in  equation  (82). 
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6.  Summary 


Conventionally,  the  Eulerian  and  Lagrangian  strain  tensors  are  derived  ei¬ 
ther  by  using  shifter  tensors  or  by  using  convected  (moving)  coordinates. 
The  definition  of  the  shifter  tensor  makes  use  of  a  scalar  product  between 
vectors  at  two  different  points  in  space  (without  first  parallel  translating 
one  vector  to  the  position  of  the  other).  When  convected  coordinates  are 
used,  vectors  and  tensors  are  associated  with  given  coordinates  in  the  con¬ 
vected  system  of  coordinates,  rather  than  being  associated  with  a  given 
point  in  the  vmderlying  space.  As  discussed  in  the  introduction,  both  of 
these  features  are  vmdesireable  when  we  need  to  understand  the  transfor¬ 
mation  properties  of  the  strain  tensors  from  an  inertial  frame  to  a  moving 
frame.  These  transformation  properties  are  also  needed  in  order  to  gener¬ 
alize  the  strain  tensor  to  Riemannian  geometry  for  appHcations  to  general 
relativity. 

I  have  provided  a  derivation  of  the  Eulerian  and  Lagrangian  strain  tensors 
for  finite  deformations  using  the  concepts  of  parallel  propagator,  the  world 
fvmction  of  J.  L.  Synge  and  the  three-point  displacement  vector  field.  This 
derivation  avoids  the  imdesireable  features  mentioned  above.  The  deriva¬ 
tion  shows  that  the  Eulerian  strain  tensor  is  a  two-point  object  that  trans¬ 
forms  as  a  scalar  under  transformation  of  material  coordinates  and  as  a  sec¬ 
ond  rank  tensor  under  transformation  of  spatial  coordinates.  The  deriva¬ 
tion  also  shows  that  the  Lagrangian  strain  tensor  behaves  as  a  scalar  under 
transformation  of  spatial  coordinates  and  as  a  second  rank  tensor  imder 
transformation  of  material  coordinates.  These  transformation  properties 
are  useful  in  understanding  how  these  strain  tensors  transform  from  one 
frame  of  reference  to  another  moving,  non-inertial  frame.  The  formulation 
presented  here  of  the  transformation  properties  of  these  tensors  is  also  use¬ 
ful  for  understanding  tire  role  of  the  reference  (xmstrained)  configuration  in 
prestressed  materials,  as  discussed  in  the  introduction. 
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